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We investigate quantum states that posses both maximum entanglement and maximum discord 
between the pertinent parties. Since entanglement (discord) is defined only for bipartite (two qubit) 
systems, we shall introduce an appropriate sum over of all bi-partitions as the associated measure. 
The ensuing definition -not new for entanglement- is thus extended here to quantum discord. Also, 
additional dimensions within the parties are considered (qudits). We also discuss nonlocality (in 
the form of maximum violation of a Bell inequality) for all multiqubit systems. The emergence of 
more nonlocal states than local ones, all of them possessing maximum entanglement, will be linked, 
surprisingly enough, to whether quantum mechanics is defined over the fields of real or complex 
numbers. 

PACS numbers: 03.67.-a; 03.67.Mn; 03.65.-w 



I. INTRODUCTION 

Entanglement, non-locality, and quantum discord are 
foundational quantum mechanical issues [T] and have re- 
ceived renewed and intense attention in the last years 
(for a very small sample, see for instance [IHl] and refer- 
ences therein). Some important distinctions between the 
three concept deserve a few words. Entanglement, non- 
locality, and quantum discord are foundational quantum 
mechanical issues [Tl and have received renewed and in- 
tense attention in the last years (for a very small sample, 
see for instance [IHl] and references therein). Some im- 
portant distinctions between the three concept deserve a 
few words. 

• While two-party entanglement is quite well under- 
stood, entanglement in a multi-partite system is an 
area of great current interest (see [51 and references 
therein). Since entanglement is a resource to be 
taken advantage of, the quest for maximally en- 
tangled states has proceeded at a rapid pace. Thus 
far, only few qubit maximally entangled states such 
as two qubit Bell states, three qubit Greenberger- 
Horne-Zeihnger (GHZ) states, and four qubit states 
[5] have been clearly identified. Of course, it would 
be of considerable interest to generate such max- 
imally entangled states as representatives of the 
ground state of a physically realizable (spin) model. 
Why? Because the strength of correlations in a 
many- body system is a reflection of the degree of 
entanglement (for pure states) [3]. We are here 
speaking of quantities such as density, magnetiza- 
tion, etc. Thus, characterization of multi- parti- 
cle entanglement and production of maximal/high 
multi-qubit entanglement is vital for the mutual en- 
richment of i) quantum information physics and ii) 
many-body condensed matter. 

• There exist certain tasks, such as device- 
independent quantum key distribution [7^ and 



quantum communication complexity problems |38j , 
which can only be carried out provided the cor- 
responding entangled states exhibit nonlocal cor- 
relations. Non-locality, as measured by the vio- 
lation of Bell's inequalities, describes the part of 
quantum correlations that cannot be reproduced by 
any classical local model. Non-locality plays a key 
role in some applications of quantum information 
theory [7] and for infinite quantum system [5], in 
which violation of Bell's inequalities can constitute 
a complementary resource and pinpoint the critical 
points associated with quantum phase transitions. 

• In turn, quantum discord [8] is defined as the dif- 
ference between two expressions of mutual infor- 
mation extended from the classical to the quantum 
realm. It aims to capture all the nonclassical cor- 
relations of a quantum system, something that en- 
tanglement fails to do. 

• In recent years, a large number of studies focus at- 
tention on strictly quantal correlations [SHTT]. and 
many different correlation measures have been pro- 
posed to detect them [TSHIS]. It is of the essence 
to point out that these correlations, and their mea- 
sures, exhibit intriguing peculiarities. Thus, viola- 
tion of Bell's inequalities implies entanglement, but 
the presence of entanglement does not necessarily 
imply Bell- violation \2A\ I25j . In turn, quantum dis- 
cord exhibits a distinctly different behavior from 
that of entanglement [26 . Remarkably enough, 
quantum discord has proved itself a useful quan- 
tity in a number of applications, such as speed-up 
in quantum computation [27_,- 

For our present purposes we specially note that quan- 
tum teleportation -one of the most intriguing features of 
entanglement- has been proposed theoretically using pro- 
tocols based on genuine multipartite entanglement, and 
experimentally realized in four and five qubits systems 
[251 m] . States with genuine multipartite entanglement 
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are defined as having strict maximal entanglement in all 
bipartitions of the system. The question of maximal en- 
tangled states and their definition was first addressed by 
Gisin PU', where the existence of only a few of them 
was first pointed out. Along lines recently pursued by 
Helwig et al. [31] , we shall focus attention here on states 
that may become maximally entangled for any number of 
parties by changing the qudit dimension. One clear mo- 
tivation for studying these states is that their maximal 
entanglement content makes them natural candidates for 
implementing original multipartite communication pro- 
tocols. 

The concept of quantum correlation for multipartite sys- 
tems is intimately related to the mathematical structure 
of quantum mechanics, as a direct consequence of the lin- 
ear character of tensor product Hilbert spaces. A ques- 
tion that will seen important below revolves around the 
fact that the Hilbert space of a system may use the field 
of real or of complex numbers. We will show instances in 
which a description provided by real coefficients is both 
useful and necessary. 



It is the aim of the present work to regard multipar- 
tite quantum correlations from three standpoints: those 
of quantum entanglement, quantum discord, and non- 
locality. 



In the Appendix we review all measures employed 
for describing quantum correlations. Some readers 
may want to peruse it before proceeding with the 
paper. Equations in the Appendix are numbered with 
a capital A preceding the corresponding number and 
are referenced throughout the text. In Section II 
we study maximally entangled states for N=3 till 8 
qubits, emphasizing the role played, in each case, by 
non-locality. Section III is devoted to obtaining maxi- 
mally entangled states for systems of qudits. States are 
compared with their corresponding qubit counterparts. 
Section IV makes use of quantum discord in order to 
find multipartite qubit system maximizing it. Section 
V compares all correlation measures by individually 
studying how they evolve as the number of parties 
increases, or globally correlating their absolute values for 
different maximal states with fixed number of parties. 
Finally, some conclusions are drawn in Section VI. 
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For simplicity, and without loss of generality, we shall 
take real coefficients (9 — in ([T])). These states max- 
imally violate the CHSH Bell inequality. Needless to 
say, only in the bipartite case entanglement implies non- 
locality and vice versa [Gisin's theorem [32 ). Let us 
recall that the family of pure states of three qubits 
|*± = (|j)±|2^=3-l-j))/y2 forms a basis, the so 
called GHZ basis or Mermin basis, and that these states 
maximally violate the Mermin inequality. In the general 
case, the most convenient basis for dealing with maxi- 
mally correlated multipartite states of N qubits will be 
given by 
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Notice that these states maximally violate the general- 
ized MABK Bell inequality. 

For N=3 qubits, the maximally entangled state is given 
by the usual GHZ state \GHZ) = (|000) + |lll))/y2. 
This instance will be the only one for which maximum 
entanglement and maximum non-locality are provided by 
the same state. For N=4 qubits, the state discovered by 
Higuchi and Sudbery [6], given by 
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with w = is the one which possesses maximum 

entanglement (9.3773 In 2) as defined by the measure of 
Eq. ( Al ). For this entanglement measure to be optimal. 



the corresponding state must possess complex expansion- 
coefficients. In addition, QD ^ for state ([3]). When 
pursuing maximum entanglement for real states only^ we 
encounter the following state 
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II. MAXIMALLY ENTANGLED STATES OF 
MULTIQUBIT SYSTEMS. THE ROLE OF 
NONLOCALITY 



When dealing with states maximizing entanglement in 
multiple qubits, it will prove convenient to employ 
the maximally correlated basis that is defined for ev- 
ery 2^-Hilbert space. Usually, the preferred basis for 
two qubit states is the so called computational one 
{|00), |01), |10), |11)}. The BeU basis of maximally corre- 
lated states for N=2 is given by 



This state-which has null discord- is highly entangled, 
but does not reach the maximal value. Instead, it stays 
very close to it (9.2017 In 2). Specifically, for the state Q 
all one qubit reduced matrices are maximally mixed, as 
in the case of the state (|3]), and two of the six associated 
density matrices for two qubits are maximally mixed. In 
the light of these results, we should require some figure 
of merit to somehow validate one of these two states or 
both of them. Let us suppose that we study the con- 
comitant MABK Bell inequality maximal violation given 
in the Appendix by Eq. (A15). In the complex case of 
four qubits, MABK^'''' = 2.17732, a 38% of the max- 
imum possible violation. On the other hand, our real 
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state of four qubits has MABK^'^ = ^6, only a 43% of 
the maximum value. This is not a big surprise, for the 
more linear combinations of maximally correlated states 
we have, the less non-locality is attained [5^. In our case, 
however, |$4^) is more nonlocal than |$4), although less 
entangled. 

In this case, the state possessing maximum correlations, 
either in the form of entanglement or of non-locality, 
is not the one with complex coefficients. Therefore, it 
makes a great difference which field of numbers we use 
for building up quantum states. 

For N=5 qubits. Brown et al. [33] proposed a state which 
is shown to possess the maximum entanglement given in 
he Appendix by (Al). This state of five qubits is of the 
form 
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with a maximum entanglement of 25 In 2. It is apparent 
that the previous state does not contain correlations go- 
ing beyond those for the bipartite case. Certainly, the 
concomitant non-locality measure is 2.1361, only 27% of 
the corresponding maximum possible value. 
For N=6 qubits, we discovered in Ref. 
a maximum entanglement given by Eq 
a state for which all bipartitions are maximally mixed 
The aforementioned state has the form 



a state with 
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Remarkably enough, Tapiador et al. |35j were able to 
algebraically simplify this state into a form that greatly 
clarifies how correlations are distributed among subsys- 
tems. This state l^'g) reads 
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It is plain from ([7| that all correlations existing in this 
maximally entangled state are encoded via maximally 
correlated 4-party and 2-party subsystems, which entails 
interesting implications regarding non-locality. Indeed, 
the maximal violation of the corresponding MABK Bell 
inequality for 6 qubits is only 2 (18% of the maximum 
value). In other words, state Q may have maximum 
entanglement, but does not violate the corresponding 
MABK inequality. For this state all concomitant re- 
duced density matrices of one, two and three qubits are 
completely mixed, yet the state is not nonlocal, j^'g) 
constitutes an example of a state which has found appli- 
cations in quantum information processing \36\ 137) , but 
since it has negligible non-locality, it may not be suitable 
for other tasks ^ i38j . 



A. Special cases: N=7 and N=8 qubits 

These two particular instances deserve special consider- 
ation for -as far as we know- they have not yet been 
investigated in detailed fashion. 

In the case of N=7 qubits, since the state obtained in 
[31] is one with no easily discernible algebraic structure, 
one wonders whether a better form for that state might 
exist, either for real or complex state-functions. A more 
detailed study is carried out here. A conjecture made 
in f31| establishes that there is no pure state of seven 
qubits with marginal density matrices for subsystems of 
one, two, or three qubits that are (all of them) com- 
pletely mixed. Contrarily to the case of the seven-qubit 
state reported in [33], we encounter other states that - 
although lacking a simple structure- all possess totally 
mixed 1- and 2-qubit marginal density matrices. Several 
highly entangled states of 7 qubits have been obtained 
(35) without reaching an upper bound. 
The aforementioned state of N=7 qubits possesses a max- 
imum entanglement of 105.7882. Certainly, its reduced 
states of two qubits are not (all) mixed, but in practice, 
they can be considered as such. This is so because, for 
reduced states of two qubits, 6 of them possess entan- 
glement 0.99205 In 2 and the remaining 15, 0.99760 In 2. 
This "imperfect" mixture surely is immaterial for practi- 
cal purposes. As far as 3-qubit marginal density matrices 
are concerned, 15 contain an entanglement of 0.97913 In 8 
and the remaining 20, 0.99541 In 8. Although some more 
compact algebraic form for the seven-qubit maximally 
entangled state could in principle be found, that is not 
actually the case. Such putative state should necessarily 
be complex -not real- to attain maximum entanglement, 
which is paramount importance for our present discus- 
sion. 

In the light of these results, we may wander if it makes 
physical sense to have states that are not completely 
mixed. Should one expect a maximally entangled state 
to be endowed with some sort of hierarchy in which in- 
creasing entanglement would be accompanied by maxi- 
mal mixedness at every one-, two-, three-stages involving 
marginal states? Also, is it mandatory to resort to states 
with complex coefficients in order to have maximum mul- 
tipartite entanglement? These are precisely the questions 
we wish to answer in the present work. 

To this end we have explored the whole space of real 
and complex states via a simulated annealing method, 
and have found similar results for both real and com- 
plex states. We have been able to detect states with 
totally mixed 1- and 2-qubit marginal density matrices, 
while increasing mixedness for 3 qubits. In the case of 
N=7 real states, we encountered a maximum entangle- 
ment of 105.1151, greater than that recorded all previ- 
ous findings. On the other hand, the complex expansion- 
coefficients instance displays a slightly greater value of 
entanglement than the real one, namely, 105.1441. While 
the former situation accrues nearly 31/35 completely 
mixed three qubit reduced density matrices, in the later 
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one that figure is increased up to 33/35. In Fig. 1, the 
coefficients of the real state are depicted. Had we let 
the optimization procedure evolve freely in the complex 
case, we would have reached the same result as in J34]. 
In other words, it is possible to obtain highly entangled 
states, either real or complex, with totally mixed 1- and 
2-qubit reduced states. To the best of our knowledge, 
no other states had been previously provided with this 
maximal entanglement values. 
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FIG. 1. (Color online) CoefHcients of the real state of A'' = 7 
qubits that majcimizes total entanglement, with fully 1- and 
2-qubit reduced states that are maximally mixed . See text 
for details. 



On the whole, although we have not been able to pro- 
vide a conclusive answer regarding entanglement for the 
N=7 qubit case, we have shed light on the characteristics 
of these states by introducing a new one that possesses 
maximally mixed reduced two qubit states. As we have 
seen, real or complex states behave differently, a situation 
that might become immaterial, provided real states suf- 
fice to fulfill eventual applications. As far as non-locality 
is concerned, all = 7 qubit instances do not violate 
the corresponding MABK inequality. Again, as in the 
N = 6 case, we find completely opposite behaviors for 
entanglement vis-a-vis non-locality. 

Now, the last case of interest, which has been quite elu- 
sive, is the one corresponding to N=8 qubits. It is undeni- 
able that as we increase more and more the dimension of 
the concomitant Hilbert space, the problem of obtaining 
an optimum state with maximum generalized entangle- 
ment becomes less and less tractable, either analytically 
or numerically. 

In any case, we have carried out extensive computations 
in order to get a suitable result. By performing a simu- 
lated annealing optimization procedure, we have reached 
a maximum for the entanglement and, in turn, obtained 
a multipartite state of eight qubits that appears to be 
an optimal one. We mention that in order to tackle the 
problem, we have only pursued real states, that is, con- 
figurations of states with real coefficients. This fact may 



limit the validity of our conclusions, but we are confident 
that, in any case, we have gotten a valid result. In ad- 
dition, we have limited the expansion-coefficients to be 
all equal in modulus, thus acquiring a simple algebraic 
structure. Needless to say, these simplifications have led 
to greater values of entanglement than arbitrary explo- 
rations with no constraints. 
Our special state of eight qubits j^g) reads 
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It exhibits interesting features. The amount of total en- 
tanglement is 362 In 2. All reduced states of one, two, 
and three qubits are maximally mixed. This fact im- 
plies that this state is perfectly suitable for i) perform- 
ing teleportation protocols involving up to three qubits 
within eight parties, as well as ii) for many other appli- 
cations. However, as expected, all four qubit reduced 
states are not completely mixed. The work by Gisin 
and Bechmann-Pasquinucci [3D] already pointed out that 
states of eight qubits cannot have all subsystems com- 
pletely mixed. However, in our case, we get close enough 
to such desideratum. 

Out of the (^) = 70 possible reduced states of four qubits, 

• 56 matrices are completely mixed (In 16) 

• 6 of them are nearly diagonal matrices with 
four equal eigenvalues (^InlG). The correspond- 
ing partitions are {1, 2, 3, 6}, {1, 2, 4, 5}, {1, 2, 7, 8}, 
{3, 4, 5, 6}, {3, 6, 7, 8}, {4, 5, 7, 8}. 

• The 8 remaining matrices are almost di- 
agonal, possessing eight equal eigenvalues 
(|lnl6). The concomitant partitions are 
{1,3,4,7},{1,3,5,8},{1,4,6,8},{1,5,6,7}, 

{2, 3, 4, 8}, {2, 3, 5, 7}, {2, 4, 6, 7}, {2, 5, 6, 8}. 

As far as non-locality is concerned, state ([8| docs not 
violate the corresponding MABK Bell inequality. In the 
light of the previous results, we may formulate the fol- 
lowing conjecture. 

Conjecture 1. No maximally entangled state greater 
that 5 qubits violates the MABK Bell inequalities. 
We shall see in forthcoming sections that non-locality for 
maximally entangled states decreases exponentially with 
the number of qubits. 

Summing up, 1) we confirm that no states of N=8 qubit 
exist with all their marginal density matrices being maxi- 
mally mixed and 2) We provided an example of real state 
that might be tailored for quantum communication pro- 
tocols or teleportation. 
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III. MAXIMALLY ENTANGLED STATES OF 
HIGHER DIMENSIONS 

Maximum qudit entanglement constitutes an extension 
of a previous study for qubits to states living in a D^- 
Hilbert space, where D stands for the dimension of each 
party. Helwig et al. [ST addressed the interesting and 
nontrivial problem of finding the conditions for the ex- 
istence of states that maximize the entanglement be- 
tween all bipartitions. These states are of the type 

\'^n,d) — 2i=o The equivalence between several 

pure state quantum secret sharing schemes and states 
with maximum multipartite (iV, D)— entanglement with 
an even number of parties is proven in [31j , an equivalence 
which indirectly implies the existence of these maximally 
entangled states for an arbitrary number of parties, based 
on known results about the existence of quantum secret 
sharing schemes. 

We will try to ascertain just what sort of states can host 
a maximum amount of entanglement between their par- 
ties. Since no general procedure has yet been provided 
for studying what conditions these maximally entangled 
qudit states should fulfill, we shall resort to numerical 
explorations that will hopefully shed some insight into 
qudits-systems that might not apply for qubits-ones. As 
dimensional-examples we have considered the cases Tif^ , 

nf^ and ^ 
The case of three qutrits {'Hf'^) is really exceptional, as 
we shall see. The Hilbert space "Hs is spanned by the 
basis {|0), |2)}. The state |*jv=3,d=3) 
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|000) - |011) - |112) + |120) - |202) + |221) (9) 



possesses maximum entanglement 3 In 3 and from inspec- 
tion, it is biseparable. The situation becomes more in- 
volved when one notices that the state 
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which is clearly nonlocal, also reaches the maximum value 
3 In 3. How is it possible that two different states (one 
being separable, and the other non-separable) may attain 
the same entanglement-based correlations? Should we 
use ( 10 1 and regard state ^ as an anomaly, or should we 
revisit instead the definition of maximum entanglement? 
Increasing individual party dimensions by one unit, we 
get a state living in the Hilbert space Tif^ (spanned 
by {|0), |1), |2), |3)}). The maximally entangled state for 
iV = 3,1? = 4 1*34) reads 
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(11) 

This dimension does not pose the puzzle we found previ- 
ously for three qutrits, in the sense that it is maximally 
entangled (3 In 4) as opposed to 
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For (111, the first two 



which has entanglement |ln4. 
elements are inseparable, while the remaining ones are 
biseparable. Thus, the role of quantum correlations is 
not as dominant as above. To be more rigorous, we 
should develop some tight Bell inequality and resort to 
the concomitant maximum violation. Unfortunately, no 
such Bell inequalities have been encountered so far. 
The case of four qutrits ("Hf"^) is the natural extension of 
the four qubit state. The state we obtain is of the form 



g ^ 10000) + |1000) + 10021) + |1021) + 
joiOO) - joilO) - 10121) - |0122)-t- 
10201) - 10202) + 10211) + 10220) 
-|1100) - |1110) + |1121) - |1122) 
-|1201) - |1202) -I- |1211) - |1220) ) 
+ ^( -|0012) + |1012) - |2001) + |2020) 

-|2102) - |2111) - |2210) + |2222) ). 

(13) 

A careful analysis shows that the above state can be 
written as a combination of tensor products, where only 
bipartite correlations appear. These loose correlations 
in the four-party case should be contrasted with the 
corresponding ones for four qubit states, where maxi- 
mum entanglement is reached for i) a state with com- 
plex expansion-coefficients (linear combination of three 
maximally correlated states) and ii) real ones (state with 
high, but not maximum entanglement embedded into a 
linear combination of six maximally correlated states). 
The differences between the cases of maximum entangle- 
ment between (iV, D = 2) and arbitrary {N, D) have to 
be taken into account in order to shed light on the prob- 
lem of quantifying and characterizing entanglement for 
multipartite systems. 

The case of five qutrits (Hf^) has turned out to be elu- 
sive. We cannot provide a simple expression for the 
concomitant state, only an approximate one. However, 
numerical evidence shows, in a 'real' quantum treat- 
ment, it is most likely that a maximum entanglement 
(25 In 3) might be reached, that is, all reduced states of 
one and two qubits are likely to be maximally mixed. 
However, this is just an approximate result. If it were 
confirmed, this would entail that in greater dimensions, 
real expansion-coefficients suffice to describe all sorts of 
states with maximum correlations. This fact has been 
confirmed for the case (N = 3,-D = 5) (5 = 3 In 5) as 
well, but the concomitant state can not be casted in sim- 
ple fashion. In view of the previous results, we formulate 
a second conjecture. 

Conjecture 2. Maximally entangled states of multiqudit 
systems {N, D) only require real expansion- coefficients. 
Also, all their reduced density matrices are completely 
mixed. 

For multiqubit states, in only a few cases all reduced den- 
sity matrices are proportional to the normalized identity 
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(that is, maximum entanglement), whereas in the case 
of qudits the instances encountered do not display this 
behavior. Furthermore, it is important to use complex 
expansion-coefRcients for qubits to get maximum entan- 
glement, but this constraint disappears in higher dimen- 
sions. These are issues of great interest that find at least 
partial elucidation in the present work. 



IV. MULTIQUBIT STATES WITH MAXIMAL 
DISCORD 



Coef ficient 


Real state 


Complex state 


Co 


0.435569236 


(0.0546370323,-0.100659299) 


Cl 


-0.186434446 


(0.0134949587,-0.188315179) 


C2 


0.151369915 


(0.465863387,0.0484337848) 


C3 


-0.0680793253 


(0.242430776,-0.0818648344) 


C4 


-0.177771681 


(0.136790716,-0.217618993) 


C5 


-0.676301307 


(-0.500174131,0.0275624382) 


C6 


-0.505730715 


(-0.442832499,0.0730590501) 


C7 


0.0567821074 


(0.0598215336,-0.379958264) 



In the case of entanglement, the most straightforward 
way of tackling quantum correlations in multipartite 
states is to introduce partitions into the system. This is 
somehow inevitable, as no definitive entanglement mea- 
sure or criterion is yet available for characterizing true 
multipartite quantum correlations of this kind. Since QD 
is defined only between pairs of qubits, it is quite natural 
to extend the same tools used for multipartite entangle- 
ment to the case of quantum discord. If not otherwise 
stated, all states are given in the computational basis 
{I000..00), I000..01), I000..10), I111..10), I111..11)}. 



A. Geometric discord 

The GQD, which is a measure introduced so as to grasp 
all properties of the usual discord measure. This is a 
computable quantity that detects (and quantifies) true 
discord. For simplicity, here we shall compute the maxi- 
mum GQD just for three qubits. 

The case of N=3 qubits is special, because the result 
we obtain is a rather simple one, but the corresponding 
state is instead of an involved nature. When approaching 
the problem of finding a state p — \<j>){4'\ of three qubits 
where GQD{pi2) + GQD {piz) +GQD{p23) is maximum, 
by definition (Cf. (AlO)), one assumes that the outcome 



for the unknown expansion-coefficients is not going to a 
simple one. Indeed, by the nature of this measure, we 
obtain via simulated annealing what appear to be ran- 
dom coefficients for the state p — \<j)){(j)\. This fact should 
not surprise anyone since the mathematical restrictions 
that are imposed when optimizing the sum of the GQD 
content for all pairs are highly nonlinear. Table [T] lists 
the real and complex coefficients -in the computational 
basis- that yield the same maximum GQD. 
The concomitant maximum GQD is equal to 5/8, but in- 
dividual discords are different! That is, GQD{pi2) — 1/8, 
GQD{pi3) = 1/4, and GQD{p23) = 1/4. Apparently, 
there is no a priori reason for this to be the case. This 
asymmetry between pairs constitutes the first discord 
feature that is different from those pertaining to the gen- 
eralized entanglement measure in multipartite systems. 
Strictly speaking, states with maximal entanglement do 
have reduced states with different entanglement values, 
but this occurs only for large number of reduced-state, 
qubit bipartitions. This sort of 'GQD-asymmetry' occurs 



TABLE I. Expansion-coefficients for the states of three qubits 
that maximize the total GQD between pairs. Columns refers 
to real and complex coefficients. See text for details. 



already in the simplest possible case of iV = 3 qubits and 
poses a conundrum. 



B. Discord measure 

Although there is an alternative measure for computing 
the discord content of a given state, we prefer the usual 
quantum discord definition QD ( A6 ) . We stress here that 



the computation of the maximum QD is a two-fold opti- 
mization procedure: first, one must find, given an arbi- 
trary state, the minimum QDs for all pairs pij , and then, 
survey all states until the maximum QD is attained. 
It would appear that since N=3 qubits is a low dimen- 
sional system, the state maximizing the concomitant QD 
should be given by an algebraically simple expression. 
As in the case of the previous GQD measure, this is 
not the case. The N=3 qubits case leads the result 
QD = 1.662026, given by the state 



1*^-°) = 0.52895|000) -F0.19492|001) -t-0.25144|010) 
-0.48224|011) + 0.38250|100) - 0.39801|101) 
-0.20213|110) + 0.20209|111). 

(14) 

The same asymmetry already found in the QGD case 
occurs also here. Also, real and complex states provide 
the same result. The non- locality of state [14] is given by 
MABKf'''' = 3.0430, a 76% of the maximum possible 
violation, with an entanglement = 2.7394 In 2 (91%, very 
close to the maximal one). Thus, a state with maximum 
QD has a relatively high amount for the several quantum 
correlations surveyed in this work. 

The N=4 case leads the resuh QD = 2.436681, with par- 
titions {1,2}, {3, 4} having QD:=l/3, and QD=0.442503 
for the remaining pairs. This symmetry also occurs in the 
case of the computation of the GQD, equal to 1.045667. 
We find the state to be of the type 
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(15) 
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FIG. 2. (Color online) QD for the family of states ( 15 1 of 



A*' = 4 qubits. Seven clear a;— values reach the maximum 
value. The a; = case is given by ( 16 1. See text for details. 



For this family of states, its QD is plotted in Fig. 2 as 
a function of x. Except for x — 0, there exists seven 
additional states that reach the maximum possible QD- 
value. The one we consider here is given in ( 15 ) with 
X = 0, that is, 



I* 
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(17) 

Notice again the tendency towards a large number of lin- 
ear combinations of maximally correlated states. Specif- 
ically, the state l^'Jgg) is particularly relevant. Again, 
complex and real coefficients lead to the same QD-value. 
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Notice that this state has all their two qubit reduced 
states of the 'X-form'. Accordingly, their QD is given in 
analytical fashion 55, • Remark that the state that we 
have found to attain maximal QD does have the same 
maximally correlated sub-states as the maximally entan- 
gled state of four qubits Needless to say, real and 
complex states reach the same maximum QD-value. Re- 
garding non-locality, MABX(|*^^)) = |%/2 (67%). 

In the case of N=5 qubits, we have not reached a "nice" 
algebraic form for the 'maximal' state. All contributions 
in the computational basis possess a nonzero weight, and 
the overall state does not seem to exhibit any symme- 
try. However, within the limits of numerical accuracy, 
all 5(5 — l)/2 — 10 pairs of qubits seem to have the 
same amount of QD, the total one being 3.642445. This 
fact implies that there ought to be a simpler form for 
the aforementioned state. The expansion-coefficients are 
shown in Fig. 3. The non-locality of the state 1^*^^) is 
« 5.0233 (63%). 

The case of N=6 qubits, on the other hand, does exhibit 
a definite symmetry. All reduced pairs have the same 
QD-value (0.350977) and the total QD is 5.264662. As 
in the case of N=4 qubits, all pairs are of the X-form and, 
thus, analytically computable. The final (real) state is of 
the form 



FIG. 3. (Color online) Coefficients of the state of = 5 
qubits maximizing QD. No clear algebraic structure can be 
drawn from these results. See text for details. 




100 120 



FIG. 4. (Color online) Similar plot as that of Fig. 3 for Af = 7 
qubits. A single state (l^'Jg 7)) significantly contributes to 
the total QD. Notice the symmetry between i and 2^ — i 
coefficients, which constitutes a clear sign of an underlying 
superposition of maximally correlated states. See text for 
details. 
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The case of N=7 qubits case is tantalizing. As in the case 
of N=5 qubits, we have not been able to find a simple 
algebraic form. Per contra, our hypothesis of -as far as 
QD is concerned- a total symmetry in the state basis is 
duly confirmed here. That is, the state of seven qubits is 
formed by a plethora (2^/2 — 64) of maximally correlated 
states. The state we obtain possesses a maximum QD of 
8.4290, which is only an approximate value. Fig. 4 shows 
the value of the expansion-coefficient for each position in 
the computational basis. The two peaks correspond to 
the state j^'^g y) , which due to unknown reasons, becomes 
differentiated from the rest. Once again, complex and 
real coefficients are equivalent when providing states with 
maximum QD. 




FIG. 5. (Color online) (a 



[N{2' 



- 1) i In 2] (upper curve). A'' = 3, and N — 5 coincide 



(1/3, horizontal line). The remaining values display a step- 
like evolution, ^jv for \GHZ]^) states (maximally nonlocal) in 
units of the same quantity [A'^(2^~^ — 1) | In 2] is also depicted 
(lower curve), (b) Value of QD^'"' vs iV. The evolution is 

N + l 

clearly exponential, (c) Normalized non- locality (to 2 2 ) for 
maximally entangled (solid line) and discordant (dot-dashed 

iV + l 

line) qubit states. The QD case saturates to fc2 2 ^ k being 
some positive constant, for an increasing number of parties, 
(d) QD'^''^ vs. S™'"'. Notice that each point corresponds 
to different states. On the whole, a monotonic increasing 
behavior is apparent. See text for details. 



V. COMPARISON BETWEEN DIFFERENT 
MEASURES 

At first sight, there is no precise way to correlate mea- 
sures S,MABK'§'^'', and QD. However, ahhough differ- 
ent in nature, they globally behave in analogous fash- 
ion. We observe that, in all three cases, the states's 
maximal measures for each quantum correlation's type 
either increase or diminish exponentially with the num- 
ber of parties. This behavior becomes apparent in Fig. 
5. Now, when relating two definite measures, special in- 
stances appear. Let us consider the case of N=4 qubits. 



It was shown that two states reached maximum "real" 
and "complex" entanglement separately. Surprisingly, 
the more entangled the state was, the lesser non-locality 
it possessed. Additionally, the state maximizing entan- 
glement is the only one in all multiqubit states here con- 
sidered that has non zero QD. In general, entanglement 
and quantum discord are more similar to each other than 
to non-locality. 

Regarding non-locality vs. entanglement, it is plain from 
Table II that maximum entanglement by no means im- 
plies maximum non-locality. Nevertheless, MABK'^'^^ 
is not that different from maximum entanglement. This 
fact is easily seen when analyzing generalized GHZ 
states, which are the ones that maximize MABK Bell 
inequalities. Measure AIABK^^^ is equal to 2^~, 
and its entanglement is SnQGHZn)) = (2^"^ - l)ln2 
{QDn{\GHZn)) = 0). Entanglement and MABK'^'''' 
for generalized GHZ states are monotonic increasing 
functions one of the other ^ Sn = (i[A/ABX]^^'=]2 - 
l)ln2. The average entanglement per partition is equal 
to In 2 in all cases. 

When computing the maximum entanglement for qubits 
systems, we obtain that 5']^'^'' > A^(2^-i-l)| ln2, equal- 
ity holding for iV = 3, 5. In Fig. 5 a) we depict S'^"'' 
and Sn{\GHZn)) in terms of [N{2'^-^ - l)iln2] vs. 
iV. What we obtain is that entanglement for maximally 
entangled states behaves in the same way, differences ap- 
pearing between clusters of states (states of 3,4, and 5 
qubits are to be compared with results for 6,7, and 8 
qubits). A kind of step-like behavior is apparent, as 
well as the decreasing tendency of Sm{\GHZm))- This 
fact constitute a strong argument concluding that non- 
locality vs. entanglement display a quite special behav- 
ior as opposed to the relation entanglement vs. discord, 
where general tendencies are more pronounced. 
Qjjmnyi Jqj. g^(.]^ number of qubits is depicted in Fig. 
5 b). This curve is perfectly fitted by an exponen- 
tial. Thus, maximum QD also evolves exponentially 
with the number of parties involved. As far as non- 
locality is concerned. Fig. 5 c) displays the evolution 
of MABK for those states with maximal discord (dot- 
dashed line) and maximum entanglement (solid line) in 

JV + l 

units of 2^^. It is plain that non-locality for maxi- 
mally entangled states diminishes exponentially, whereas 
in the case of quantum discord, a sort of proportionality 
is reached {MABK{QD'^'''') cx 2^ for TV > 7). 
When comparing maximum quantities corresponding to 
different states, we obtain a result such as the one de- 
picted in Fig. 5 d). That is, the QD for those states 
maximizing QD monotonically increases with their en- 
tanglement counterpart. 



VI. CONCLUSIONS 

We have explored the tripod entanglement-discord-non- 
locality for multiqubit systems. We have revisited the 
known cases of maximally entangled states of iV = 4- 



9 



iV 


MABK{S^a^) 


MABKiQD^a.) 


Si^QDrnax ) 


3 


4 - 100% 


3.0430 - 76% 


2.7394 In 2 - 91% 


4 


2.1773 - 38% 


|%/2- 67% 


8 In 2 - 85% 


5 


2.1381 - 27% 


^ 5.0233 - 63% 


17.7864 In 2 - 71% 


6 


2 - 18% 


5^2 - 62.5% 


41.2542 In 2 - 63% 


7 


2 - 13% 


^ 7.2524 - 45% 


« 97.9457 in 2 - 64% 


8 


2 - 9% 







TABLE II. Several types of correlation measure for states 
with increasing number of qubits. The first column displays 
de maximum violation of the generalized MABK Bell inequal- 
ity for those states that maximize the generalized measure of 
entanglement, as well as the concomitant percentage with re- 

N + l 

sped to the maximum possible violation (2 2 for each N). 
The second row displays similar results for those states with 
maximum QD. Finally, in the last column, the entanglement 
for states maximizing QD is computed. Notice that some re- 
sults are analytic. The computation of QD for states with 
maximum entanglement is not depicted, as it vanishes except 
for N=4 {QD{Sma^) = 0.7548 - 31%). See text for details. 



, 5-, and 6-qubits, exploring their properties as far as 
non-locality is concerned. Specifically, the cases of 7 
and 8 qubits have been studied in some detail, and al- 
lowed for original results in both cases. It has become 
clear that, some similarities notwithstanding, entangle- 
ment and non-locality MABK'^^^ behave differently in 
one or the other direction. That is, the more entangled 
the states are, the less nonlocal they become. Recall 
that, however, we only focus on extreme cases, and not 
on arbitrary states. The aforementioned tendency has 
been confirmed for an increasing number of qubits. This 
is the case because entanglement and non-locality consti- 
tute, after all, different resources. 

Quantum discord, which has been recently proved to be a 
valuable resource, has been extended here to multiqubit 
systems by introducing a measure over all two-qubit par- 
titions. Employing the corresponding maximum value, 
we have discovered states of "maximum discord" for up 
to 8 qubits. These states possess very interesting fea- 
tures, including unexplained asymmetries among their 
parties. Furthermore, we strongly believe that these 
maximal QD— states may find interesting applications in 
quantum information processing. 

In the light of recent research, we have also explored 
entanglement in multiqudit systems, and obtained the 
concomitant maximally entangled states for different 
{N, £))— configurations. Also, the way of defining quan- 
tum mechanics over real or complex numbers is subject to 
inspection here. We find that real expansion-coefficients 
in a given basis suffice to explain maximum entangle- 
ment. Recall that this has not been the case for multi- 
qubit states. 

Finally, the comparison of all types of measure leads us 
to conclude that entanglement and discord are both pos- 
itively correlated, whereas non-locality decreases with ei- 



ther maximum entanglement or discord. We conclude, 
thus, that maximum entanglement and maximum dis- 
cord display some similarities, whereas this scenario is 
blurred when considering non-locality. 



Appendix A: On correlations measures 

Research on the properties and applications of multipar- 
tite entanglement measures has attracted considerable 
attention in recent years ^.33, 39 48 . One of the first 
useful entanglement measures for A^-qubit pure states 
10) to be proposed was the one introduced by Meyer 
and Wallach [35]. It was later pointed out by Brennen 
[40] that the measure advanced by Meyer and Wallach 
is equivalent to the average of all the single-qubit linear 
entropies, that is, the average entanglement of each qubit 
of the system with the remaining (A'^— l)-qubits. 

Another way of characterizing the global amount of 
entanglement exhibited by an A^-qubit state is provided 
by the sum of the (bi-partite) entanglement measures 
associated with the 2^"-'^ — 1 possible bi-partitions of 
the A^-qubits system |33j . This particular number takes 
into account that the marginal density matrices describ- 
ing the kth party, after tracing out the rest, are equiv- 
alent to those oi N — k parties because of the relation 
(^) = (jv^^)- In essence, these entanglement measures 
are given by the degree of mixedness of the marginal den- 
sity matrices associated with each bi-partition. In our 
case, we shall use the von Neumann entropy 



SvN = ^ -Tr[pi\npi], 



(Al) 



where the sum is performed over all 2^~^ — 1 different 
bipartitions. 

A considerable amount of research has been devoted to 
unveil the mathematical structures underlying entangle- 
ment, in particular concerning those states which possess 
maximum entanglement, as given by some appropriate 
measure. The survey of (pure) states maximizing mea- 



sure ( Al ) has been the subject of intense research, initi- 
ated by the work of Brown et al. [33] . These genuinely 
entangled states achieving the maximum possible entan- 
glement according to the ( Al ) have recently investigated 



with regards to their usefulness for perfect quantum tele- 
portation, superdense coding, and quantum secret shar- 
ing. Indeed, highly entangled multipartite states gener- 
ate intense interest for quantum information processing 
and one-way universal quantum computing ^49j . They 
are essential for several quantum error codes and com- 
munication protocols [501, since they are robust against 
decoherence. The aforementioned research for maximally 
entangled states will be extended here to N—7 and N—8 
qubits in Section III. 

Quantum discord is also a valuable resource for the imple- 
mentation of non-classical information processing proto- 
cols [2^ [5m54) . In the light of these developments, it be- 
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comes imperative to characterize multipartite states not 
only through their maximum amount of entanglement 
but also through the maximization of their correspond- 
ing total quantum discord, establishing the links within 
the tripod entanglement-discord-nonlocality. 
In addition to entanglement, quantum discord [HI [22] con- 
stitutes a rather recent information-theoretical measure 
of the "non-classicality" of bipartite correlations as given 
by the discrepancy between the quantum counterparts of 
two classically equivalent expressions for the mutual in- 
formation. Quantum discord corresponds to a new facet 
of the "quantumness" that arises even for non-entangled 
states. More precisely, quantum discord is defined as the 
difference between two ways of expressing (quantum me- 
chanically) such an important entropic quantifier. Let p 
represent a state of a bipartite quantum system consist- 
ing of two subsystems A and B. If S{p) stands for the 
von Neumann entropy of matrix p and pA amd ps are 
the reduced ( "marginal" ) density matrices describing the 
two subsystems, the quantum mutual information (QMI) 



Mq reads W 



Mqip)^SipA) + S{pB)^S{p). 



(A2) 



This quantity is to be compared to another quantity 
Mq{p), expressed using conditional entropies, that clas- 
sically coincides with the mutual information. To define 
Mq{p) we need first to consider the notion of conditional 
entropy. If a complete projective measurement 11^ is 
performed on B and (i) pi stands for Tvab P and (ii) 
Pyi||n^ for [n^ pUf /pi], then the conditional entropy be- 
comes 



out numerically, unless the two qubit states belong to the 
class of the so called X-states, where QD is analytic [55) . 
Accordingly, we shall introduce here a generalized QD 
measure for multiqubit states in the form 



Ptl 



(A6) 



where the sum takes place only over all N{N — l)/2 re- 
duced two qubit states pi, since QD is only defined in 
that case. 



Since the evaluation of QD ( A6 ) involves an optimiza- 



tion procedure and analytical results are known only in 
a few cases, an interesting alternative was advanced in 
[35] by introducing a geometric measure of quantum dis- 
cord (GQD). Let X be a generic zero QD state. The GQD 
measure is then given by Hilbert-Schmidt norm 



GQDip)^Mm^[\\p~x\ 



(A7) 



where the minimum is taken over the set of zero-discord 
states X- Given the general form of an arbitrary two- 
qubits state in the Bloch representation 

4(0 = I (g) I + Y.l=l ^uCTu ®1 + X]«=l Vu^ ®(Ti + 

(A8) 

with Xu = Tr{p{au®T)), y„ = Tr(p{I(^au)), and Tuv = 
Tr{p{au ® cTv)), it is found in Ref. ^Hj that a necessary 
and sufficient criterion for witnessing non-zero quantum 
discord is given by the rank of the correlation matrix 



5(A|{nf}) = ^p,5(p^linf), 



and Mq{p) adopts the appearance 



(A3) 



Mq{p) 



{nf} 



S{pA)~S{A\{Tlf}). (A4) 



Now, if we minimize over all possible 11^ the difference 
Mq{p) — Mg(/9)|jjB-|. we obtain the QD, that quantifies 
non-classical correlations in a quantum system includ- 
ing those not captured by entanglement. The most gen- 
eral parameterization of the corresponding local measure- 
ments that can be implemented on one qubit (let us call it 

B) is of the form {11% = /^(g) |0') (0'|, = 7^(8) |1')(1'|}- 
More specifically we have 



/ 1 2/1 2/2 2/3 

xi Til Ti2 Ti3 

X2 T21 T22 T23 

\ T31 T32 



(A9) 



that is, a state p of the form ( A8 1 exhibits finite quantum 



discord iff the matrix (A9| has a rank greater that two. 



It is seen that the geometric measure ( A7 1 is of the form 



gqd{p) = Um^ + \\t 



J?, 



(AlO) 



|0') ^ cosa'lO) -he^^'sina']!) 
|1') ^ e""'' sin a' 1 0) - cos a' 1 1) 



(A5) 



which is obviously a unitary transformation -rotation in 
the Bloch sphere defined by angles (a',/3')- for the B 
basis {|0),|1)} in the range a' G [0, tt] and 13' <E [0,27r). 
The previous computation of the QD has to be carried 



where ||x|p — J2u^u^ Kimx is the maximum eigen- 
value of the matrix (a;i, 0:2, a;3)*(a;i, X2, X3) + TT'^ and 
R — 1/Tr(p^). The study of the properties of the gen- 
eralized GQD will be here carried out for N ~3 a nd N=4 
qubits, its form being analogous to that of (A6). 



Entanglement, however, is not the only quantum corre- 
lation that can be considered in multipartite quantum 
systems (qubits in our case). Since the formalization by 
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Werner [53] of the modern concept of quantum entangle- 
ment it has become clear that there exist entangled states 
that comply with all Bell inequalities (BI). This entails 
that non-locality, associated to Bl-violation, constitutes 
a non-classicality manifestation exhibited only by just a 
subset of the full set of states endowed with quantum 
correlations. 

In some cases, however, entangled states are useful to 
solve a problem if and only if they violate a Bell in- 
equality |38j . Moreover, there are important instances 
of non-classical information tasks that are based di- 
rectly upon non-locality, with no explicit reference to 
the quantum mechanical formalism or to the associated 
concept of entanglement [7]. Most of our knowledge on 
Bell inequalities and their quantum mechanical viola- 
tion is based on the CHSH inequality [S7]. The scenario 
with two dichotomic observables per party, is the sim- 
plest one [58 endowed with a nontrivial Bell inequal- 
ity for the bipartite case (with binary inputs and out- 
comes). Quantum mechanically, these observables reduce 
to Aj(Bj) = aj(bj) • (T, where aj(bj) are unit vectors in 
M'^ and a — (ctx, Cy, cr^) the Pauli matrices. Violation of 
CHSH inequality requires the expectation value of the op- 
erator BcHSH = Ai(g)Bi-f Ai(g)B2-f A2(8)Bi-A2«)B2 
to be greater than two. It is indeed an inequality that 
poses an upper limit (2-\/2, the Tsirelson-bound [W), to 
quantum mechanical correlations between distant events. 
A relation exists here with the discussion and experimen- 
tal determination of whether local variable model (LVM) 
variables are required for, or even compatible with, the 
representation of experimental results. A proper way to 
measure non-locality (for the two qubit state p) uses the 
quantity 



omax 
^CHSH 



maxTr{pBcHSH) 

a;,bi 



(All) 



Similarly, non-locality in the three qubit case is given by 
the violation of the Mermin's inequality |60| . This in- 
equality was conceived originally in order to detect gen- 
uine three-party quantum correlations, impossible to re- 
produce via LVMs. The Mermin inequality states that 
Tr{pBMerrnin) < 2, whcrc BMermin IS the Mermin oper- 
ator 



B„ 



B, 



aib2b3 ' 



B, 



-^6162031 (A12) 



where B^vw = u- cr(X'v-cr(g)w-o' with a — (a^, (^yi<yz) 
being the usual Pauli matrices, and aj and bj unit vectors 
in 'M? . Notice that the Mermin inequality is maximally 
violated by Greenberger-Horne-Zeilinger (GHZ) states. 
As in the bipartite case, we shall define the following 
quantity 



Mermin" 



max Tr{pBM e 







(A13) 



as a measure of the non-locality of the state p. While in 
the bipartite case the CHSH inequality is the strongest 



possible one, no bond equivalent to Tsirelson's is avail- 
able for arbitrary dimensions. 

Now, in the case of multiqubit systems, one must instead 
use a generalization of the CHSH inequality to N qubits. 
This is done in natural fashion by considering an exten- 
sion of the CHSH or Mermim inequality to the multipar- 
tite case. The first Bell inequality (BI) for four qubits was 
derived by Mermin, Ardehali, Belinskii, and Klyshko j60j . 
One deals with four parties with two dichotomic out- 
comes each, the BI being maximum for the generalized 
GHZ state (|0000) + |llll))/\/2. The Mcrmin-Ardehali- 
Belinskii-Klyshko (MABK) inequalities are of such na- 
ture that they constitute extensions of older inequalities, 
with the requirement that generalized GHZ states must 
maximally violate them. To concoct an extension to the 
multipartite case, we shall introduce a recursive relation 
that will allow for more parties. This is easily done by 
considering the operator 



Bn+, (X [{B, + B[) (E)Bn + (Si - B[) (E) B^], (A14) 

with _Bjv being the Bell operator for N parties and Bi = 
V • (T, with a = {(Tx, ay, Gz) and v a real unit vector. The 
prime on the operator denotes the same expression but 
with all vectors exchanged. The concomitant maximum 
value 



MABKI^'''' = max Tr{pBN) 

ai,bi 



(A15) 



will serve as a measure for the non-locality content of a 
given state p of N qubits if aj a nd bj are unit vectors in 
M'^. The non-locality measure (A15) will be maximized 



by generalized GHZ states, being the corresponding 
maximum value. 

The MABK inequalities are not the only existing Bell 
inequalities for N qubits [BT, but they constitute a sim- 
ple generalization of the CHSH one to the multipartite 
case. Accordingly, it will suffice to use these particular 
inequalities to illustrate the basic results of the present 
work, as far as non-locality is concerned. 
Regardless of the quantity one shall compute or of the 
state one aims at finding, some sort of optimization pro- 
cedure must b e ca rried out. In the case of maximum 
entanglement (All, the 2^ coefficients of the multipar- 



tite pure state constitute the variables to deal with (twice 
if the state is of complex instead of real nature). For the 
quantum discord QD, a minimization takes place for the 
two parameters in every N{N — l)/2 reduced two qubit 
state, whereas non-locality MABK'j^^^ requires an ex- 
ploration among their corresponding unit vectors defin- 
ing the observers' settings. In any case, we have per- 
formed a two-fold search employing i) an amoeba opti- 
mization procedure, where the optimal value is obtained 
at the risk of falling into a local minimum and ii) the 
well known simulated annealing approach [62! . The ad- 
vantage of this computational 'duplicity' is that we can 
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be confident regarding the final results reached. Indeed, knowledge partial support under project FIS201 1-23526 
the second recipe contains a mechanism that allows for 
local searches that eventually can escape local optima. 
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